NATIONAL  ADVISORY  COMMITTEE 
FOR  AERONAUTICS 

TECHNICAL  NOTE  2372 


RECTANGULAR-WIND-TUNNEL  BLOCKING  CORRECTIONS  USING 
THE  VELOCITY -RATIO  METHOD 
By  Rudolph  W.  Hensel 

Southern  California  Cooperative  Wind  Tunnel 
California  Institute  of  Technology 


Washington 
June  1951 


Best  Available  Copy  I  200008)6  )90  I 

DISTRIBUTION  STATEMENT  A  QUAim"  4 

Approved  for  Public  Release 
.  Distribution  Unlimited  /tAMoo-/! 


NATIONAL  ADVISORY  COMMITTEE  FOR  AERONAUTICS 


TECHNICAL  NOTE  2372 


RECTANGULAR-WIND-TUNNEL  BLOCKING  CORRECTIONS  USING 
THE  VELOCITY-RATIO  METHOD 
By  Rudolph  W.  Hensel 


SUMMARY 


In  this  report  calculations  of  the  ratios  of  the  velocity  increments 
at  test  bodies  to  those  at  the  tunnel  walls  caused  by  the  solid  blocking 
of  these  bodies"  within  the  walls  of  a  closed  rectangular  wind  tunnel  are 
presented.  The  bodies  treated  include  two-dimensional  airfoils;  small 
bodies  of  revolution;  straight ,  untapered,  finite-span  wings  of  varying 
span;  and  swept,  untapered,  finite-span  wings  of  varying  span.  It  is 
shown  that,  after  wake  blocking  effects  have  been  removed,  the  present 
method  furnishes  semiempirical  blocking  corrections  for  most  wind-tunnel 
models  and  their  components.  Results  are  presented  for  all  the  cases 
mentioned.  The  test-section  proportions  of  the  Southern  California 
Cooperative  Wind  Tunnel  at  the  California  Institute  of  Technology  (viz., 
ratio  of  height  to  width  equal  to  1^)[2)  are  used  in  calculations. 


INTRODUCTION 


The  velocity— ratio  method  of  obtaining  blocking  corrections  in 
high-speed,  subsonic  wind  tunnels  was  first  solved  hy  Gothert  for  the 
cases  of  a  body  of  revolution  and  of  a  finite-span  wing  of  span-to- 
diameter  ratios  of  0.23  and  0.3  in  a  closed  circular  tunnel  (reference  1). 
This  work  was  later  extended  in  an  unpublished  report  to  the  cases  of  a 
wing  having  a  span— to— diameter  ratio  of  0.75  nnd  of  a  wing  spanning  a 
closed  circular  wind  tunnel. 

In  the  present  report,  the  methods  of  reference  2  are  used  to  extend 
the  previous  results  for  straight  wings  of  varying  span  in  a  closed  rec¬ 
tangular  tunnel  whose  height-to-width  ratio  is  l/f2.  For  swept  wings,  a 
slightly  different  approach  involving  the  use  of  line  doublets  has  been 
utilized. 

This  work  was  done  by  the  California  Institute  of  Technology  and 
has  been  made  available  to  the  NACA  for  publication  because  of  its 
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SYMBOLS 
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c 

d 

g 

H 

I 

M 

Mu 


m 

n 

q 

q’ 

r 

R 

s 

s'  =  s/W 


distance  of  point  source  and  sink  from  doublet  origin 
center,  or  axis  of  test  section 
wing  chord 

maximum  diameter  of  a  body  of  revolution 
x-coordinate  of  a  source  segment  , 

tunnel  height 

length  of  a  body  of  revolution 
Mach  number,  corrected  for  blocking 

tvuinel  calibration  Mach  number,  uncorrected  for  blocking 
(may  include  strut  calibration) 

integer,  indicating  image  number  in  y-direction 

'  I 

integer,  indicating  image  number  in  z-direction 

line-source  strength,  square  feet  per  second 

point-source  strength,  cubic  feet  per  second 

distance  from  source  or  doublet  element  to  point  at  which 
velocity  increment  shall  be  obtained;  +  y^  + 

remainder 

semispan  of  model  wing  (measured  in  y-direction) 


b 


t 


maximiim  wing  thickness 


11  total  axial  velocity  increment  in  test  section  due  to  all 

images  except  primaiy  one  (mo'del) 

u'  total  axial  velocity  increment  in  test  section  due  to  all 

images^  including  primary  one  (model) 

axial  velocity  increment  in  test  section  due  to  a  single 
image 

U  tunnel  axial  velocity 

velocity  in  x-direction  due  to  a  point  doublet 
W  tunnel  width 

w  radial  coordinate;  =  y^  + 


X  axial  coordinate 

Xi  =  X 

y  lateral  coordinate 


z  ^  vertical  coordinate 
Zp  =  z  -  nH 

9  =  sin“^  S. 

r 

line-doublet  strength,  cubic  feet  per  second 
M''  point-doublet  strength,  feet^  per  second 

A  angle  of  sweep  of  a  given  wing  at  any  Mach  number  M^ 

Aq  equival^t  angle  of  sweep  at  M^  =  0 

y  ^  ^  .Stokes  stream  function,  cubic  feet  per  second 


u 
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Subscripts : 

A,  C,  and  particular  points  in  test  section  or  on  test-section  „ 

so  forth  walls 

For  further  explanation  of  the  symbols  see  figures  1  to  U. 

DETERMINATION  OF  INTERFERENCE  VELOCITIES 
Two-Dimensional  Wing 


A  two-dimensional  wing  may  be  represented  by  a  chordwise  distri¬ 
bution  of  infinite  line  sources  and  sinks.  The  axial  velocity  increment 
produced  by  any  single  infinite-line-source  image  (fig.  l)  is 


^  -  q 

sin  9 

qg 

(1) 

1  2Tir 

2nT^ 

at  a  point  A  in  the  center  of  the 

basic 

tunnel , 

-  n^H^  +  g^.  Thus, 

=  qg/  1  \ 

^A  2n^j^2jj2  ^  g2y 


(2) 


Omitting  the  central  source  and  svtmming  for  the  remainder  gives 


(3) 


Ik 


Similarly,  the  axial  velocity  at  the  wall  of  the  tunnel,  point  B, 
including  both  the  central  source  and  all  the  images,  is 
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For  norraal  chord  sizes, 
and  hence  may  be 

with  the  result 


(5) 


Thom  has  shown  in  reference  2  that  qg  for  a  single  line  source  may  be 
replaced  by  ^qg,  the  distribution  of  sources  and  sinks  representing 
the  wing  section.  Therefore,  since  the  ]^qg  terms  would  also  cancel 
each  other,  the  result  obtained  for  a  single  line  source  in  equation  (5) 
is  identical  to  that  for  a  complete  wing. 


g  is  small  compared  with  the  tunnel  dimensions 
neglected  in  comparison  with  n^  or  ^n  -  , 


Bocfy  of  Revolution 

A  body  of  revolution  may  be  represented  a  distribution  of  point 
sources  and  sinks  along  the  tunnel  center  line.  The  axial  velocity 
increment  due  to  a  single  image  point  source  (fig.  2)  is 


ui  = 


q' 

Unr^ 


sin  9 


q'g 

Unr3 


(6) 


at  point  A  in  the  center  of  the  basic  tunnel. 
Then,  substituting  into  equation  (6), 


u-. 


q*g 

Unw3 


2  IT 


3/2 


(7) 


Again  omitting  the  central  source  and  summing  for  the  remainder  with  the 

2 

assumption  that  g^/w^  is  negligible  compared  with  2^  and  m^  gives 

W^- 
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Similarly,  it  can  be  shown  that  the  wall  velocities 


and  u*Q  are 


q'g  ^ 

00  00 
uEZ 

n=l  m=l 

1 

3/2  » 

' 

n=l 

1 

-V 

3/2 

L  .  if  ]£  t  ,.2 

(n  -  if  i£ 

> 

(m=0) 

2/  W2j 

j 


and 


-  q'g< 

00  00 

n=l  m=l 

1 

3/2  ^ 

.2^ 
m=l 

1 

3/2 

c  ,  31 

l+nW^ 

2  1)2  /  1\2 

L  if 

/ 

i 

L"  ?  ("  ■  2)J 

(„.0) 

Lr  ■  5)J 

(10) 


The  latter  two  equations  include  the  effect  of  the  primaiy  source,  which 
must  be  omitted  in  the  calculation  of  u^^. 


I 


Straight,  Untapered,  Finite-Span  Wing 

A  finite-span  wing  may  be  represented  by  a  distribution  of  finite- 
length  sources  and  sinks.  The  axial  velocity  increment  produced  by  a 
single  source  element  is  (fig.  3),  for  this  case. 


du^ 


q  <^72. 


sin  0 


qg  dyi 


(11) 


For  the  point  A  in  the  center  of  the  basic  tunnel,  the  general  expres¬ 
sion  for  the  square  of  the  distance  from  the  source  element  is 

^A^  =  +  n^H^  +  ^mW  +  Integration  across  the  image  span  gives 

the  total  contribution  at  A  of  one  image;  namely. 


NACA  TN  2372 


7 


p  /  - - 

''  [g2  -  * 


+  ^mW  + 


Performance  of  the  integration  leads  to  the  result 


mW  +  s 


mW  -  s 


l|7i(n^H^  +g^)j[g2 

+  +  (mW  +  s  ,2]V3 

[g2  +n2H2  +  (mW-s)2]' 


which  is  the  same  as  equation  I6,  reference  2.  Making  equation  (13) 

e2 

nondimensional  and  again  neglecting  the  —  terms  give 


qg  /  1  ' 

,  »rT.r2 1  ^  tt2 


S 

^  “  W 


■A  Itnw  l  2I£  r  2  nfi  '  f  2 


("  - 1) 


As  before  for  the  case  of  the  body  of  revolution,  the  total  velocity 
increment  is 


q2_2_  +  22_+22_ 

n=l  m=l  m=l  n=l  A 


However,  it  will  be  noted  that  the  single  summation  for  which  n  =  0 

leads  only  to  the  indeterminate  quantity  This  may  be  evaluated  by 

the  application  of  L'  Hospital's  Rule,  finally  giving  for  a  single 
source  line 
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Similarly,  the  total  axial  velocity  increments  at  points  B  and  C  on 
the  walls  of  the  tunnel  are 
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9 


As  before,  the  primary  source  line  which  corresponds  to  the  wing  is  not 
included  in  calculating  the  velocity  increments  at  any  point  in  the  wing. 
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Swept,  Untapered,  Finite-Span  Wing 

In  the  preceding  treatment  for  unswept  bodies  the  analyses  were  * 

carried  out  for  single  line  and  point  sources  as  a  simplification  for 
bodies  which  could  be  represented  by  a  combination  of  line  and  point 
sources  and  sinks.  Thus,  for  a  given  value  of  g,  all  the  image 
sources  are  in  the  same  plane  and  a  constant  distance  from  the  plane 
of  A,  B,  and  C,  namely,  g.  Obviously,  this  condition  no  longer  holds 
for  swept  bodies  in  which  g  would  be  a  function  both  of  span  and 
angle  of  sweep,  and  hence  the  simplified  treatment  is  no  longer  appli¬ 
cable.  It  is  then  necessary  to  proceed  to  a  more  general  representation 
of  the  flow  field  involving  the  use  of  both  sources  and  sinks.  One  of 
the  simplest  of  such  combinations  is  the  doublet  and  that  is  what  is 
used. 


The  Stokes  stream  ftinction  for  a  point  doublet  is  given  by 
(reference  3) 


i  = 


(20) 


2aq  * 

where  [j,'  =  - - .  The  point-doublet  strength  p.'  remains  finite  as  a. 

hn  ’ 

the  distance  of  the  point  source  and  sink  from  the  doublet  origin, 
approaches  zero  and  q'  goes  to  infinity.  The  velocity  in  the  x-direction 
caused  by  the  doublet  element,  whose  axis  is  parallel  to  the  x-axis 
(fig.  1;),  is  given  by 


Then,  for  a  swept-wing  element,  with  doublet-element  axes  parallel  to 
the  x-axis. 


dy 


1 


duj  =  [i 


(22) 


NACA  TN  2372 


11 


The  square  of  the  distance  from  the  doublet  element  to  A  is 

^A^  ~  ^1^  n^H^  +  ^mW  +  Considering  Aq  positive  for  sweepback, 

then  for  the  right  wing  with  respect  to  the  model,  —  =  tan  Aq,  and 

-X-j_  ^1 

for  the  left  wing,  -  =  tan  Aq.  Having  eliminated  x-]_,  the  expres- 

^1 

Sion  for  the  contribution  of  a  single-image  doublet  is 
%  ■  “  £ 


y^^tan^AQ  +  n%^  +  ^mW  +  y^^ 


3/2 


3  n^H^'  +  ^mW  + 


y^^tan^AQ  +  n^H^  +  ^mW  +  y^j ^ 


5/2 


dy- 


(23) 


It  should  be  noted  that  a  change  of  sign  between  the  right  and  left 
halves  of  a  wing  makes  no  difference  mathematically.  Therefore,  the 
swept  wing  acts  in  exactly  the  same  manner  as  a  wing  yawed  at  an 
angle  Aq.  Performing  the  integration  and  letting 


.  =  n2  Hi 

^  W2 

b^'  =  2m 


c  =  tan'^AQ  +  1 


dl'  =  W'c  -  (bA*)^ 


S  * 

_  s 

W 

Xa.' 

II 

+ 

bA-s’ 

+  c(s' )2 

V 

11 

1 

CO 

+  c(s' 
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In  following  a  similar  procecure  to  obtain  and  it  is  found 

that  the  constants  are  the  same  as  for  u.  except  that  the  subscripts 

^  1 

are  changed  to  B  and  C,  respectively,  and  n  — >  n  -  —  at  B  and 

1  ^ 
m — >  m  -  ^  at  C.  The  summations  are  the  same  as  for  the  straight 

wing;  namely, 
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Support  Struts 

The  method  previously  used  to  calculate  the  velocity  ratios  for 
aerodynamic  bodies  may  also  be  applied  to  support  struts.  However,  in 
order  to  avoid  infinite  velocity  increments  at  the  junctions  of  struts 
and  model,  it  is  necessary  to  consider  the  support  system  in  the  tunnel 
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as  an  integral  part  of  the  model.  With  this  viewpoint,  velocity  ratios 
could  also  be  calculated  for  struts  in  the  same  manner  as  for  wings,  in 
which  just  the  images  are  summed  in  order  to  calculate  u^.  The  span 

of  the  struts  would  be  constant  at  2s/W  =  1.0  in  a  reoriented  tunnel 
for  which  h/W  =  \[2. 

Using  points  A  and  C,  where  C  would  now  be  at  the  top  of  the  reori¬ 
ented  tunnel,  the  velocity  ratio  %/n'c  for  half  wings  or  struts  would 

be  identical  to  those  for  complete  wings  or  double  struts,  both  com¬ 
pletely  spanning  the  tunnel.  There  is  one  simple  case  for  which  the 
result  is  immediately  known,  namely,  a  single,  centrally  mounted, 
unswept  strut  for  which  the  two-dimensional  result  of  l/3  applies  (see 
section  Two-Dimensional  Wing).  In  general,  if  other  than  a  single-strut 
support  system  were  used,  the  velocity  ratio  would  be  a  f\mction  of  the 
strut  spacing  used  as  well  as  of  the  angle  of  sweep.  The  problem  would 
be  further  complicated  by  the  presence  of  a  rear  strut,  which  is  fre¬ 
quently  the  case. 

Because  of  the  additional  complexity  involved  and  the  expected 
difficulty  in  separating  the  total  wall  velocity  increment  into  the 
separate  effects  due  to  the  solid  blocking  of  the  support  system,  model 
wing,  and  model  fuselage  and  to  the  wake  blocking  of  each  component,  all 
of  which  may  have  different  velocity  ratios,  no  general  solution  of  the 
•support-strut  case  has  been  presented.  It  would  be  sin^jler  and  probably 
more  accurate  to  perform  a  complete  calibration  of  the  area  in  which  a 
model  would  normally  be  mounted,  with  the  struts  installed.  During  this 
calibration,  the  wall  pressures  at  B  and  C  could  also  be  obtained,  thus 
giving  base  values  which  include  both  the  solid  and  wake  blocking  and 
interference  effects  of  the  model  support  system. 


NUMERICAL  CALCULATIONS  AND  RESULTS 


The  methods  used  in  summing  the  doubly  infinite  series  are 
explained  in  detail  in  appendix  A.  Briefly,  calculations  are  made  for 
each  image  up  to  a  finite  nmber  n  =  m  =  n-j^.  The  remainders  are 

obtained  by  direct  integration  from  n^  to  infinity  after  making  certain 

simplifying  assTunptions.  Unfortunately,  the  series  convergence  is  not 
very  rapid  and  it  is  necessary  to  take  n]^  as  high  as  seven  in  most 
cases.  Fortunately,  however,  there  is  a  negligible  difference  between 
the  remainder  terms  for  the  swept  and  straight  wings,  since  the  effect 
of  sweep  rapidly  diminishes  as  the  distance  between  images  and  tunnel 
increases.  Hence  the  very  diffic\ilt  problem  of  attempting  to  integrate 
the  complicated  remainders  for  the  swept  wings  is  avoided.  In  all 
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cases,  the  dimensional  ratio  of  the  Southern  California  Cooperative 
Wind  Tunnel,  namely,  H/W  =  l/\f2,  has  been  used. 

Figures  5  and  6  show  the  variation  with  2s/W  of  the  velocity 
ratios  and  Uj^/u'q,  respectively,  for  sweep  angles  of  0°,  30°, 

il3°,  and  60°.  It  should  be  noted  that  each  set  of  curves  begins  from  a 
common  value  at  2s/W  =  0,  which  corresponds  to  the  body-of-revolution 
case,  and  that  the  effect  of  sweep  is  to  increase  the  ratios  over  the 
values  for  Aq  =  0.  The  effect  of  the  side  walls  upon 

%/^'C  "the  straight  wings  as  2s/W  approaches  unity  is  also  of 

some  interest.  In  the  former  case,  as  the  wing  approaches  the  side 
walls  their  effect  upon  ^a/^'b  continues  until  the  wing  just  touches 

the  walls,  at  which  point  the  ratio  discontinuously  drops  from  the 
three-dimensional  result  which  includes  the  effect  of  the  side  walls  to 
the  purely  two-dimensional  value  of  l/3  in  which  the  side-wall  effect 
has  been  eliminated.  As  the  wing  tip  nears  the  side  wall  at  point  C 
(2s/W— ^1)  it  also  causes  u'q  to  approach  infinity  while  uj^  remains 
finite,  so  that  uWu'q  approaches  zero. 


It  should  be  noted  that,  although  the  velocity  ratios  for  swept 
wings  are  greater  than  those  for  straight  wings,  the  absolute  magnitudes 
of  u^  and  u'g  decrease  with  increasing  sweep  for  a  given  span.  The 

value  of  u'q  is  also  less  for  Aq  /  0  than  for  Aq  =  0,  but  behaves 

irregularly  beyond  2s/W  =0.7.  These  three  blocking  velocity  incre- 

ments,  made  nondimensional  by  multiplication  with  are  plotted 

in  figures  7,  8,  and  9  as  a  function  of  wing  span.  The  numerical  values 
given  for  these  wings  are  based  upon  line  doublets.  The  corresponding 

J I  ttW*^  iiTiW^ 

values  of  uj^  ,  and  u'q  for  a  body  of  revolution. 


which  is  represented  by  a  point  doublet,  are  16.2,  56.8,  and  29.k, 
respectively.  From  figure  9  it  can  be  seen  that  u'q  becomes  negative 


for  small  angles  of  sweep  as  the  wing  tip  approaches  the  wall 
(2s/W->  1.0). 


In  reference  2,  Thom  has  used  an  arbitrary  line  or  point  source  of 
strength  q  displaced  a  distance  g  from  the  origin,  giving  the 
multiple  qg  which  appears  in  his  equations.  Thom  considered  two 
elementary  shapes  in  two  and  three  dimensions,  namely,  a  Rankine  oval 
and  ovoid,  and  the  shapes  corresponding  to  a  line  and  point  source  of 


strength  +q  at  a  distance  g  from  the  origin,  a  line  and 
strength  ^  q  at  the  origin,  and  a  line  and  point  sink  of 


point  sink  of 
strength  ^  q 
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at  a  distance  -g  from  the  origin.  For  the  Rankine  oval  and  ovoid, 
which  are  the  bluntest  bodies  obtainable  by  this  method,  yqg  is  2qg, 

)l  0 

which  must  be  a  maximum.  For  the  latter  cases,  )  qg  =  ^  qg.  Thom  has 

calculated  g/c  and  ^  as  a  function  of  thickness  ratio  t/c  for 

Ut  Q 

both  two-dimensional  shapes  and  g/l  and  — tt  as  a  function  of  d/l 

Ud'^ 

for  both  body-of-revolution  shapes.  This  procedure  may  be  repeated  for 
any  arbitrary  source-and-sink  distribution  to  facilitate  a  more  correct 
calculation  of  %/U,  u'g/U,  and  u'q/U  by  Thom's  method  for  a  given 

shape  of  wing  or  fuselage. 

Although  the  numerical  values  of  figures  7,  8,  and  9  are  based 
upon  line  doublets,  they  have  been  divided  by  and  hence  are  equiva¬ 
lent  to  the  results  for  a  line  source  divided  by  qg,  since  p  =  2qg. 

Therefore  Thom's  method,  wherein  ^qg  for  a  specific  shape  replaces 
qg  or  p  when  using  these  figures,  is  still  directly  applicable. 

Using  the  values  of  ^  and  g/c  given  in  figure  1;  of  reference  2 

for  a  Rankine  oval,  one  can  easily  calculate  the  incompressible-flow 
values  of  u^^/U,  u'g/U,  and  u'q/U.  For  example,  for  a  straight  wing 

of  this  section  having  unit  chord,  thickness  ratio  t/c  =  0.1,  and 
span  2s/W  =  0.8  in  a  tunnel  of  width  W  =  12  feet,  =  0.00075. 


APPLICATION  OF  METHOD 


In  the  preceding  sections,  ratios  of  the  velocities  at  the  model 
to  those  at  reference  points  on  the  wall  have  been  calculated  for 
certain  representative  bodies.  During  a  given  test,  then,  the  velocities 
at  the  walls  must  be  measured  in  the  model  plane  ABC,  and  so  forth,  and 
compared  with  the  reference  values  obtained  without  the  model  installed. 
The  difference  of  the  two  is  not  immediately  u'g  or  u'q,  however. 

Included  in  the  measurement  are  the  effects  of  the  struts ,  if  the  model 
is  held  in  this  manner  and  if  the  base  wall  pressures  do  not  include 
the  strut  effect,  and  also  an  amount  due  to  wake  blocking.  It  will  be 
assumed  that  the  tunnel  has  been  calibrated  with  the  struts  installed. 

The  wake  blocking  presents  a  more  complicated  problem,  but  its  magnitude 
is  fortunately  small  compared  with  the  solid  blocking  for  most  models. 

The  methods  of  separating  the  wake  and  solid  blocking  effects  and 
obtaining  the  final  corrected  dynamic  pressure  and  Mach  number  M  from 
the  measured  wall  pressures  have  been  outlined  in  reference  J4.  The 
procedure  given  therein  for  the  separation  of  wake  and  solid  blocking 
may  be  somewhat  modified  for  straight  wings  and  fuselages  by  the  use  of 
the  wall  pressures  rather  than  use  of  the  theoretical  wake  blocking  as 
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given  by  Thom  (reference  2).  The  latter  has  represented  the  wake  of  a 
body  of  revolution  or  a  wing  by  a  point  or  line  source^  respectively^ 
at  the  body  with  a  corresponding  point  or  line  sink  downstream  at 
infinity.  If  the  sources  are  assumed  to  be  in  the  model  plane,  their 
axial  contribution  in  that  plane  is  zero.  This  constitutes  a  good 
assumption  for  a  body  of  revolution  and  a  straight  wing.  Thus  for  these 
cases  the  only  wake  effects  at  the  model  plane  are  due  to  the  sinks  at 
infinity,  which  create  a  uniform  disturbance  in  the  flow  field.  There¬ 
fore,  the  velocity  ratios  due  to  wake  blocking  are  unity  for  straight 
wings  and  three-dimensional  bodies.  The  absolute  magnitude  of  the  wake 
blocking  effect  may  be  obtained  by  the  fact  that  the  solid  blocking 
effects  at  the  walls  disappear  downstream  of  the  model  and  the  remaining 
velocity  increment  approaches  the  total  disturbance  due  to  the  point  or 
line  sources  and  sinks  representing  the  wake.  The  wall  wake  increment 
at  the  model  plane  is  then  just  one-half  of  the  total  far  downstream 
because  of  the  zero  contribution  o:^  the  sources  at  the  model  plane. 
Therefore,  the  wall  wake  increment  may  be  determined  and  subtracted  from 
the  total  velocity  increments  at  B  and  C,  enabling  the  solid  blocking 
corrections  to  be  evaluated. 


The  experimental  procedure  outlined  above  is  not  strictly  appli¬ 
cable  to  the  case  of  swept  wings  because  they  extend  downstream  from 
the  reference  plane  ABC,  which  only  passes  through  the  apex  of  the 
wings.  Therefore,  the  velocity  ratios  due.  to  the  wake  will  no  longer 
be  unity  because  of  the  finite  contribution  made  by  the  wing  source 
elements  to  the  wall  velocities  at  B  and  C.  However,  for  wings  of  small 
sweep  and  of  small  span  compared  with  the  tunnel  width  combined  with 
fuselages,  for  which  the  previous  argument  still  holds,  it  may  be 
expected  that  the  velocity  ratios  due  to  wake  blocking  in  the  apex 
plane  ABC  are  sufficiently  close  to  unity  to  permit  use  of  this  procedure 
with  sufficient  accuracy. 


After  the  strut  and  wake  effects  upon  the  wall  velocity  have  been 
eliminated,  the  velocity  ratios  given  in  this  report  may  be  applied. 
However,  the  factors  for  wing  and  fuselage  (body  of  revolution)  are  in 
general  dissimilar,  and  it  is  necessary  to  estimate  the  relative  magni¬ 
tude  of  the  correction  due  to  each.  This  may  be  done  from  the  model 


solid  blocking  results  of  reference  2  at  the  model  for  fuselages  and 
straight  wings  or  by  comparing  the  results  of  figure  7  with  the  body- 

of -revolution  result  of  u^  =  17.0,  knowing  ^qg  for  both  wing 


and  fuselage;  the  estimated  relative  magnitudes  at  the  wall  may  be 
obtained  by  application  of  the  respective  velocity  ratios.  As  mentioned 


previously,  sweep  decreases  the  velocity  increments  at  all  points  in  the 
flow  field.  Figure  10  shows  the  variation  of  the  ratio  of  Uj^  with 


sweep  to  Vl^  for  Aq  =  0  as  a  function  of  Aq  for  various  spans. 
This  figure  will  facilitate  estimation  of  the  relative  magnitudes  of 
the  fuselage  and  wing  corrections  when  ^  0°. 
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Now  that  the  measured  values  of  u'/U  due  to  solid  blocking  at 
points  B  and  C  on  the  walls  have  been  split  up  into  the  separate  contri¬ 
butions  of  the  wing  and  fuselage,  u/U  at  point  A  in  the  model  for 
each  is  known.  However,  the  increase  in  velocity  is  generally  of  lesser 
interest  than  the  corresponding  increase  of  Mach  number  and  dynamic 
pressure.  For  small  corrections,  the  equations  given  by  Thom  (refer¬ 
ence  2)  and  also  in  reference  i;,  revised  to  the  present  notation,  are 


The  ratio  q/q^  is  now  the  ratio  of  the  dynamic  pressure  corrected  for 
solid  blocking  to  the  uncorrected,  or  calibration  dynamic  pressure. 

Using  the  above  corrections  and  a  calculated  value  of  the  wake 
blocking  for  a  given  model,  it  is  possible  to  operate  at  a  desired 
corrected  Mach  number  at  the  model  by  setting  a  predetermined  wall  Mach 
number  at  point  B  or  C .  This  procedure  is  particularly  desirable  in 
tunnels  where  rapid  determination  of  the  test  results  is  required,  such 
as  at  the  CWT. 

The  velocity-ratio  method  of  obtaining  blocking  corrections  has 
certain  advantages  over  previously  used  methods.  First  of  all,  being 
semiempirical  in  nature,  it  promises  greater  accuracy,  particularly 
since  the  velocity  ratios  are  usually  less  than  unity.  Hence  the  cor¬ 
rected  Mach  numbers  and  dynamic  pressures  are  less  than  the  measured 
ones,  whereas  by  other  methods  they  are  higher,  being  the  result  of 
increases  of  lower  observed  values  by  theoretical  correction  factors. 
Procedures  which  raise  observed  Mach  numbers  to  higher  corrected  values 
by  means  of  theoretical  formulas  are  always  inaccurate  and  dangerous  at 
subsonic  speeds  closely  approaching  M  =  1.  However,  correcting  down 
from  a  high  observed  Mach  nmber  to  a  lower  corrected  value  as  is  done 
with  the  wall-pressure  procedure  is  much  safer  and  more  accurate. 
Secondly,  as  Gothert  has  pointed  out  in  reference  1,  the  Prandtl-Glauert 
rule  no  longer  applies  at  the  body  in  the  supercritical  speed  range  but 
is  again  applicable  some  distance  from  the  model.  In  fact,  Gothert 
states  that  the  velocity-ratio  method  may  be  applied  with  good  accuracy 
up  to  the  speed  at  which  the  shock  waves  approach  the  tunnel  wall, 
namely,  the  choking  speed.  Most  blocking  theories  use  the  Prandtl- 
Glauert  rule  and,  of  course,  break  down  as  the  critical  Mach  number  of 


NACA  TN  2372 


21 


the  model  is  exceeded.  Furthermore,  although  the  present  analysis  has 
been  carried  out  at  zero  lift,  Gbthert  has  pointed  out  in  reference  1 
that  the  procedure  is  applicable  to  a  first  approximation  to  cases  with 
lift  by  averaging  the  pressures  at  the  two  walls  opposite  the  lifting 
surfaces.  This  would  normally  be  at  points  B  and  E. 


The  derivations  of  the  preceding  sections  have  been  carried  out 
using  line  and  point  sources  and  doublets.  These  unit  singularities 
may  be  replaced  by  combinations  of  sources  and  sinks  properly  dis¬ 
tributed  to  represent  the  desired  shape.  The  factors  qg  and  p, 
determine  the  body  shape  for  straight  wings  and  bodies  of  revolution 
and  for  swept  wings,  respectively.  Replacement  of  these  factors  by. 
^qg  eind  ^p,,  which  are  chosen  to  simulate  the  desired  shape,  leads 
to  the  actual  velocity  increments  in  the  incompressible-flow  field  for 
this  particular  body.  However,  it  is  seen  that  these'  shape  factors  are 
common  multiples  in  all  equations  for  each  case;  therefore  the  velocity 
ratios  are  independent  of  shape. 


It  is  also  shown  in  reference  1  that  a  given  size  of  compressible- 
flow  field  corresponds  to  a  fictitious  incompressible  one  in  which 
the  axial  (x)  dimension  is  unchanged  and  the  lateral  and  vertical 


^c  \ll  -  \  and 


(y  and  z)  dimensions  are  shrunk  by  an  amount  y^^ 

Zi  =  Zq  \|l  -  ,  where  the  subscripts  c  and  i  refer  to  com¬ 

pressible  and  incompressible,  respectively.  Therefore,  for  straight 


wings  and  bodies  of  revolution,  ^qg  and  ^p,  are  functions  of  Mach 
number,  but  the  velocity  ratios  in  the  plane  of  the  origin  (x  =  O) 
remain  constant.  This  also  applies  to  swept  wings  as  far  as  airfoil 
shape  is  concerned  but  the  application  of  this  rule  results  in  a  change 
of  plan  form.  In  the  fictitious  incompressible  plane,  the  span 
decreases  while  the  axial  distances  remain  constant,  causing  the  angle 
of  sweep  to  increase.  The  physical  sweep  angle  A  of  a  given  wing  at 
a  given  value  of  M^  in  the  compressible  plane  then  corresponds  to  a 
value  of  sweep  for  M^  =  0,  namely,  Aq,  which  is  required  by  the 
present  method.  Figure  11  shows  the  variation  of  Aq/A  with  M^^  for 

different  values  of  A.  Interpolation  in  figures  5  and  6  will  give  the 
corrv.ct  velocity  ratio  at  any  Mach  number  for  a  given  wing. 


There  are  some  limitations  of  this  blocking-correction  procedure 
which  must  also  be  considered.  It  will  be  noted  that  in  each  case  the 
axial  distance  between  the  sources  and  the  origin  or  between  the  sources 
and  sinks  approaches  zero.  Therefore,  it  has  tacitly  been  assumed  that 
the  body  length  or  chord  must  be  small  compared  with  the  dimensions  of 
the  tunnel.  Should  such  not  be  the  case,  a  small  correction  must  be 
applied  to  the  velocity  ratio.  The  magnitude  and  nature  of  this  cor¬ 
rection  for  a  circular  wind  tunnel  are  given  in  reference  1.  For  the 
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small  models  normally  encountered  in  high-subsonic-speed  testing,  this 
correction  may  usually  be  neglected.  The  derivations  have  also  been 
carried  out  for  constant-chord  wings,  while  most  models  have  tapered 
wings.  The  effect  of  taper  may  be  taken  into  account  by  assuming  a 
spanwise  variation  of  g,  but  this  would  involve  much  additional  alge¬ 
braic  complexity,  particularly  for  swept  wings,  which  would  multiply 
the  amount  of  work  in  the  series  siommations.  Furthermore,  only  the 
velocity  ratios  at  the  centers  of  the  wings  have  been  calculated,  while 
the  mean  along  the  span  should  be  used  for  any  given  wing.  For  straight 
wings,  the  additional  calculations  could  easily  be  carried  out  using 
equation  (19).  For  swept  wings,  however,  each  position  along  the  span 
is  in  a  different  axial  plane  and  would  require  an  additional  set  of 
equations  of  greater  complexity.  The  errors  involved  in  neglecting  the 
effects  of  taper  and  of  the  spanwise  variation  of  velocity  ratio  tend 
to  compensate  each  other. 

In  conclusion,  it  should  be  mentioned  that  these  results  may 
easily  be  applied  to  another  type  of  testing  frequently  employed  at  the 
Cooperative  Wind  Tunnel.  This  is  the  semispan-reflection-plane  method 
wherein  the  model  is  mounted  on  the  tunnel  floor  directly  without  any 
strut  support  system  and  with  the  wing  extending  vertically  from  the 
floor.  It  is  to  be  noted  that  the. procedure  simply  doubles  the  effec¬ 
tive  test-section  area  but  does  not  alter  the  dimensional  ratio 
H/W  =  l/\f^.  Therefore,  the  preceding  results  are  directly  applicable 
to  a  reoriented  test  section  in  which  A  is  the  center  line  of  the  tvuinel 
floor  on  which  the  model  is  mounted,  B  and  E  are  in  the  floor  corners, 
and  C  is  then  in  the  center  line  of  the  tunnel  ceiling. 


Southern  California  Cooperative  Wind  Tunnel 
California  Institute  of  Technology 
Pasadena,  Calif.,  May  21,  19ij.8- 
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APPENDIX  A 

CALCULATION  OF  REMII'TOERS 


As  mentioned  previously,  the  convergence  of  the  series  derived  in 
the  text  of  this  paper  is,  in  general,  quite  slow.  Thus  for  practical 
purposes  it  is  necessary  to  carry  the  numerical  calculations  for  the 
contributions  due  to  individual  images  to  some  reasonable  values  of  n 
and  m,  say,  n  =  m  =  n^^.  Furthermore,  if  it  can  be  assmed  that  n]_ 

is  large,  certain  simplifying  assumptions  can  be  applied  to  the  series 
equations,  enabling  a  direct  integration  in  both  the  m-  and  n-dlhections 

to  be  performed  from  n2_  +  ^  to  infinity.  The  value  n^^  +  i,  inci¬ 
dentally,  comes  about  from  the  fact  that  the  individual  series  summation 
for  n^  extends  to  the  outer  boundary  of  that  particular  image  and 
hence,  since  the  point  n  =  m  =  n^  is  in  the  center  of  the  image,  the 
outer  boundary,  at  which  the  integration  must  begin,  is  at  n^  +  1. 

Two  cases  are  considered,  namely,  that  of  the  body  of  revolution 
in  the  center  of  the  tunnel,  and  the  straight,  untapered,  finite-span 
wing.  It  is  readily  apparent  that  the  effect  of  sweep  rapidly  becomes 
negligible  as  the  distance  from  the  primary  image  becomes  large  compared 
with  the  dimensions  of  a  swept  wing.  Since  n^  is  assumed  large,  the 

remainder  terms  for  swept  wings  may  be  considered  to  be  identical  to 
those  for  straight  wings.  The  very  rapid  convergence  of  the  swept  to 
the  straight  case  is  immediately  apparent  in  the  calciilations  for  indi¬ 
vidual  images  for  which  m  and  n  are  less  than  n^.  A  value  of  seven 
was  used  for  n^^  in  all  numerical  calculations. 

In  all  cases  the  remainders  have  been  broken  into  separate  parts 
in  a  similar  manner.  Since  n^^  is  assumed  to  be  large,  it  follows 

that  the  remainders  are  independent  of  the  particular  location  in  the 
primary  image  (i.e..  A,  B,  G,  etc.)  under  consideration.  Hence,  it  is 
sufficient  to  calculate  the  results  for  only  one  quadrant  and  simply 
multiply  by  four.  In  the  actual  calculations ,  an  integration  was 

carried  out  from  n]_  +  |  to  infinity  in  both  the  n-  and  m-directions 

and  the  result  miltiplied  by  four.  Then  the  image  strips  for  m  =  0 

were  integrated  from  np  +  |  to  infinity  in  the  n-direction  and  multi- 

plied  by  two.  This  was  repeated  for  the  individual  strips  for  m  =  1 
through  n^  in  the  n-direction  and  the  results  for  each  were  multiplied 
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by  four.  The  same  procedure  was  carried  out  for  the  individual  image 
strips  in  the  m-direction.  The  above  may  be  readily  seen  in  the 
following  analyses. 


,  Body  of  Revolution 

The  velocities  at  A,  B,  and  C  due  to  the  solid  blocking  effect  of 
a  body  of  revolution  in  the  center  of  a  rectangular  wind  tunnel  are 
given  by  equations  (8),  (9),  and  (10),  respectively.  Using  the  pro¬ 
cedure  outlined  in  the  preceding  paragraphs,  the  results  are 
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i;nW^ 

q'g 


n=l  ni=l 


The  remainder  terms  are  the  double-integral  term  and  the  two  braced 


quantities, 


,  and  are  identical  for  A,  B,  and  C 


under  the  assumptions  used.  Performance  of  the  indicated  integration 
results  in 
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Therefore,  the  remainder  is  a  fiinction  only  of  n^  and  the  dimensional 

ratio  of  the  tunnel  H/W  and  may  readily  be  calculated  for  any  par¬ 
ticular  case. 


Untapered,  Finite-Span  Wing 

As  mentioned  previously,  the  effect  of  sweep  upon  the  remainder 
as  n-j_  becomes  large  is  negligible.  Therefore,  the  simpler  equations 

for  Aq  =  0  may  be  integrated  and  applied  to  the  swept-wing  case  as 

well.  The  velocities  at  A,  B,  and  C  due  to  the  solid  blocking  effect 
of  a  straight,  finite-span  wing  are  given  by  equations  (I6),  (17),  and 
(18),  respectively.  Again  applying  the  assumption  that  n2_  is  large, 
these  equations  may  be  transformed  into 
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The  remainder  terms  again  are  the  double-integral  term  and  the  two 


braced  quantities, 


and 


R 


n 


R 


and  are  identical  for  A,  B, 


m 


and  C  under  the  assumptions  used.  Performance  of  the  indicated  inte¬ 
gration  leads  to  the  result 


Since  the  remainder  for  the  wing  case  is  a  function  only  of  H/W, 

and  s/W,  it  may  be  calculated  for  any  particular  case. 
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Figure  4.-  Swept  finite-span  wing 


revolution 


Figure  5*-  Variation  of  u. /u'p,  with  span. 


H/W  =  l/\}2 


Angle  of  sweep,  Aq,  deg 

Figure  10.-  Variation  of  ratio  of  u^  with  sweep  to  without  sweep 

with  angle  of  sweep. 


Mach  number, 


Figure  11.-  Variation  of  Aq/A  with  Mach  number. 


tan  Aq  = 


tan  A 
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semiempirical  blocking  corrections  for  most  wind-  semiempirical  blocking  corrections  for  most  wind- 

tunnel  models  and  their  components.  The  test-  tunnel  models  and  their  components.  The  test- 


